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Regularization of the Covariant Derivative on Curved Space 

by Finite Matrices 

Masanori Hanada*) 
Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

\& ' In a previous paper [ M. Hanada, H. Kawai and Y. Kimura, Prog. Theor. Phys. 114 

(2005), 1295 ] it is shown that a covariant derivative on any n-dimensional Riemannian 
manifold can be expressed in terms of a set of n matrices, and a new interpretation of IIB 
matrix model, in which the diffeomorphism, the local Lorentz symmetry and their higher spin 
j-j \ analogues are embedded in the unitary symmetry, is proposed. In this article we investigate 

several coset manifolds in this formulation and show that on these backgrounds, it is possible 
to carry out calculations at the level of finite matrices by using the properties of the Lie 
algebras. We show how the local fields and the symmetries are embedded as components 
of matrices and how to extract the physical degrees of freedom satisfying the constraint 
proposed in the previous paper. 



! 81. Introduction 

vo . 

Although it is conjectured that string theory provides the unification of fun- 
damental interactions, its present formulation based on perturbation theory is not 
satisfactory. In order to examine whether it actually describes our four-dimensional 
world, a non-perturbative and background independent formulation is needed. Ma- 
trix models represent a promising approach to studying the nonperturbative dy- 
namics of string theory. For a critical string, they are basically obtained through 
j^H. dimensional reduction of the ten-dimensional U (N) N = 1 supersymmetric Yang- 

J-h ; Mills theory. 2 )- 4 ) 

IIB matrix model 3 ) is obtained through dimensional reduction to a point, and 
the action is given by 

^. 

. S = -\Tr(\[A a ,A b ][A\A b ] + U 1 a [A a ^]\, (1-1) 



where V is a ten-dimensional Majorana-Weyl spinor, and A a and ip are N x N 
Hermitian matrices. The indices a and b are contracted by the flat metric. This 
action possesses an SO(10) global Lorentz symmetry and U(N) symmetry. One 
problem with this model is that it is unclear how curved spaces are described and 
how the fundamental principle of general relativity is realized in it. 

It may seem that the space of dynamical variables becomes very small after 
the dimensional reduction, but in fact this is not the case if N is infinitely large. 5 ) 
Indeed, in Ref. 1) a new interpretation, in which IIB matrix model contains gravity 
in a background independent manner, is given. The argument presented there is 
as follows. The matrix variables A a and ij) a act on the Hilbert space V = C N as 
endomorphisms, i.e. linear maps from V to itself. Because V is infinite dimensional in 

* ) hanaQgauge . scphy s . kyot o-u . ac . j p 
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the large- AT limit, we can give various different interpretations to it. If we assume V is 
a space consisting of an n-component complex scalar field, i.e., V = {(p 1 : R 10 — ► C ra }, 
instead of V = C N , an endomorphism T is a bilocal field K l ^(x,y), which can be 
formally regarded as the set of differential operators of arbitrary rank with n x n 
matrix coefficients: 6 ) 

n ,. 
3=1 J 

= jfaixWi?) + c$(x)d^(x) + ^f(x)d fl d^(x) + • • • . (1-2) 

In particular, we can take the covariant derivative as a special value of A a , 

A a = i(d a -ia a (x))eEnd(V). (1-3) 

In this sense, the ten-dimensional U (n) J\f = 1 super Yang-Mills theory can be 
embedded in the matrix model, and local gauge symmetry is realized as a part of 
the U (N) symmetry of the matrix model, 

5A a = i[\,A a ], (1-4) 

where A is a matrix- valued function on M 10 , which is a 0-th order differential operator 
in Endiy). With this interpretation of V, however, the model does not contain 
gravity. In order to embed gravity, we regard A a to be a covariant derivative acting 
on a curved space. Then, the diffeomorphism and the local Lorentz symmetries 
become parts of the U(N) symmetry of the matrix model. In this interpretation, 
any curved space corresponds to a certain matrix configuration, and the path integral 
includes the summation of all the curved spaces. 

The procedure for describing the covariant derivative with a set of matrices is 
the following. 1 ^ Let M be an n-dimensional Riemannian manifold with a fixed spin 
structure and M = Ujf/j be its open covering. On each patch U{, the covariant 
derivative is expressed as 

V^ = e a m ^(d m -u> m hc ^O bc ), (1-5) 

where e a m and to m bc are the vielbein and the spin connection, respectively. Here, 
Ob c is the Lorentz generator that acts on Lorentz indices. The index [i] is the label 
of the patch. In the overlapping region UiHUj, Vk J and Vk JJ are related as 

vW = Rj>(t tJ (x))V^, (1-6) 

where tij : Ui n Uj — > G = Spin(n) is the transition function and R a b (tij(x)) is the 
vector representation of Uj(x). 

Let us consider the principal G bundle on M associated with the spin structure, 
and denote it by E pr i n (M). It is constructed from the set of Ui x G by identifying 

(x\i],9[i\) with (x {j] , g {j] ): 



x [i]= x [jh 9\i] = Uj(x)g\j}. 



(1-7) 
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We take V = C°° (E pr i n (M)) , which is the space of smooth functions on E pr i n (M). 
We assume that covariant derivatives act on the space V; that is, O ab generates an 
infinitesimal left action, 

ie ab (0 ah f®) (x, g) = /« (x, (l + ie ab M ab ) g\ - /« (x, 5 ) , (1-8) 

where M ab is the matrix of the fundamental representation. Then, we can construct 
endomorphisms from a covariant derivative as follows: 

V[1) = «(.) 6 ^ 1 )V?. (1-9) 

Here, R^ b (g) is the vector representation of G.*) Each component of Vr a ) is globally 
defined on E pr i n (M) (i.e. vK = V/^ in an overlapping region) and is indeed an 
endomorphism on V . The index (a) merely labels n endomorphisms. 

This method is valid in any number of dimensions, and we can express the co- 
variant derivative on any n-dimensional Riemannian manifold in terms of n matrices. 

Based on the procedure described above, the new interpretation of IIB matrix 
model 1 ) mentioned above becomes feasible. If the matrices A a are sufficiently close 
to the covariant derivatives iV( a ) on one of the manifolds M, it is natural to regard 
each A a as acting on C°° (E pr i n (M)) and to expand A a about iV( a ):**) 

K = iV( ) + a {a) (x,g) + ia (a) {b) (x, g)V {b) + ia {a) bc (x, g)O bc 

+i 2 a (a) W(x,g)V (b) V {c) + i 2 a {a) ^ d ( x ,g)V [b) O cd + -- - . (Ml) 

In this expansion, local fields appear as coefficients. For example, a^ b \x, g) con- 
tains the fluctuations of the vielbein. Coefficients of higher-derivative terms corre- 
spond to fields of higher spin. In this sense, this part of the space of large- iV matrices 
describes the dynamics around the background spacetime M. The diffeomorphism 
and the local Lorentz symmetry are realized as parts of the U(N) symmetry of the 
matrix model, which are generated by A = §{A^ a \ V( a )} and A = i\ ab O a b***^ Note 
that we do not fix M and all Riemannian manifolds with all possible spin structures 

*' R a b and R( a ) b are the same quantity. However, we formally distinguish them, because the a 
and (a) obey different transformation laws. Specifically, a is transformed by the action of G, while 
(a) is not. 

**' Strictly speaking, because A a is Hermitian, we should introduce the anticommutator { , } in 
Eq. itTTB : 

A a = iV (a) + a (a) (x,g) + l{a (a) w (x,g),V (b) } + ^{a {a) bc {x, g),O bc } + ■■■ . (1-10) 

***' Similar arguments have been given in the case of the spectral action. 7 - 1, 8 - ) In this case, V is 
taken to represent the spinor fields, and only bispinor operators, such as the Dirac operator, are 
considered. This is sufficient to include the symmetries of general relativity, although we cannot 
treat the action 11-11 in this formalism. By comparing it with our formalism, it may possible to 
clarify how the physical degrees of freedom are embedded in matrices. 
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are included in the path integral. In Ref. 9), this formalism is extended to the su- 
permatrix model, so that the local supersymmetry is included in the superunitary 
symmetry 

Now let us consider the equation of motion. 1 ) Variation with respect to A a gives 
the equation of motion 

[A a ,[A a ,A b ]] = 0. (1-12) 
Here we simply set ip = 0. If we impose the ansatz 

K = iV( a ), 

Eq. 1)1-12(1 becomes 



= 0. (144) 
Note that Eq. (|144|) is equivalent to 

= [v<Mv a ,v 6 ]] 

= [V a ,R ab cd O cd ] 

= (v a R ab cd ) O cd - R b c V c , (145) 

where R a b cd is the Riemann tensor and R b c = R a b ac is the Ricci tensor. We have 
assumed that V a is torsionless. Eq. Q145|) holds if and only if 

V a R ab cd = 0, R ab = 0. (146) 

The first equation here follows from the second by the Bianchi identity \/^ a R bc de ^ = 0. 
Therefore, the covariant derivative on a Ricci-flat spacetime is a classical solution. 
Recently, it is also shown that the equation of motion of the massless higher spin 
gauge fields about flat space can be derived. u > 

In general, the method presented in Ref. 1) can be applied only in the large-iV 
limit, and we must introduce some regularization. In order to perform a Monte-Carlo 
simulation, or to calculate a divergent quantity like the free energy, an explicit reg- 
ularization using finite- N matrices is necessary. This is also useful for the purpose 
of extracting information concerning the topology from the matrices. The exten- 
sity of the eigenvalue distribution does not reveal the topology of the corresponding 
manifold; indeed, any manifold entails an infinitely large momentum. In order to 
determine the topology, we must investigate finer structures, e.g. commutation re- 
lations and the degeneracy of eigenvalues. On generic backgrounds, it is difficult 
to write the explicit forms of regularizations in terms of finite- TV" matrices, but for 
certain manifolds with large symmetries, it is possible. In this paper, we give a few 
examples: spheres, real and complex projective spaces, tori and flat spaces. 

The organization of this paper is as follows. In £JH we discuss the heat kernel 
regularization**^. Although this method is not so useful in actual calculations, it can 

*' The same action has been studied in a different context. 10 ' 

**' This section is partiaily based on collaboration with H. Kawai and Y. Kimura. 
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be applied to any background, and it makes clear the meaning of the regularizations 
presented in £j31 I n 32]we introduce the explicit regularization procedures for spheres, 
real and complex projective spaces, tori and flat spaces. On these backgrounds, the 
corresponding E pr i n have the structures of Lie groups, and by using their algebraic 
properties, explicit calculations become tractable. We investigate the case of S 2 in 
detail and show how the local fields and symmetries are embedded as components 
of matrices. Section is devoted to conclusions. In Appendix El we give formulae 
used in £J2 In Appendix [Bl we discuss the classical solutions of modifications of IIB 
matrix model with mass or cubic terms. 



§2. Heat kernel regularization 



In the previous section, we introduced the general procedure to express the co- 
variant derivative in terms of a set of matrices. The Hilbert space V = C°°(E pr i n (M)) 
is a space of smooth functions on E pr i n (M). In order to regularize the trace in a 
general covariant way, the heat kernel regularization 12 ) is useful. 

Let A pr i n be the Laplacian on E pr i n (M), which is defined by 

A prm = Y,V {a) 2 + K 2 Y,°ab 2 , (2d) 

(a) a,b 

where k is an arbitrary constant with dimensions of mass that specifies the mass 
scale in the direction of Spin(n). As k becomes larger the damping factor on the 
higher spin fields becomes more stringent. The heat kernel regularization is defined 
by 

TrX = J dg J ed d x(x, g\X\x, g) 

j dg j ed d x(x,g\e tA P" n/2 Xe tA P" n/2 \x,g), (2-2) 



Tr t X = K d - D 



where f is a parameter with dimensions is (mass) 2 , dg is the Haar measure of 
Spin(n), e = dete m a is the determinant of the vielbein, and D = d ^ d ^ \ s the 
dimension of E pr i n (M). Eventually we take the limit t — * 0. Using the heat kernel 
regularization, we can cut off the higher frequency modes of C oc {E pr i n {M)). In terms 
of the fields on M, the modes with large momentum or large spin*) are suppressed. 
Expanding the matrix variables in local fields as Eq. (jlTlj) and using the heat kernel 
regularization, we can express physical quantities in terms of the local fields. As an 
example, let us evaluate the bosonic part of the action. If we ignore the local fields 
other than the graviton, then using a formula given in Appendix^J it becomes 

Sbosonic = ~4~2 Tr * ([ iV (a)> iV (fe)] 2 ) 
1 „ /„ \ 2 



Ag 2 



2 Tr t {R(a)(b) Cd O c d 



*^ Note that the momentum along the direction of Spin(n) correspond to the spin. 
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l K d - D -H- D l 2 - X 



jdgj ed d xR abcd R abcd + O^/ 2 " 2 ). (2-3) 



4g 2 (4tt) d / 2 

If we include the higher spin fields, their kinetic terms and nonlinear interactions 
involving derivatives appear. For example, consider the case in which the spin 3 field 
a a mn (x) is excited. For simplicity, let us assume that the background spacetime is 
flat. Then we have 

Sbosonic = -J^2 Tr t ( i9 (a) + \ { a (a) (C)W , 3(c) 9{d) } , »3(6) + \ { a (b) ^ , } 
1 K ^ t -D/2-3 



4<? 2 4(4tt) d / 2 



dg I ed d x(a a cd {d c a a ^)a a c ' d '{d c a h ^') 

-a a cd {d c a a e S)a b c ' d \d c a a e 'f) 



x (<W<W%' + 14 permutations) + 0{r D/2 ~ 2 ). (2-4) 

The kinetic term of a a bc (x) appears at 0{t~ D l 2 ~ 2 ). In general, the action becomes 
far more complicated, as an infinite number of higher spin fields appear and interact 
with each other. Couplings with the curvature of the background spacetime also 
exist. 

Some remarks are in order here. First, although the heat kernel regularization 
preserves the general covariance, it breaks the higher spin gauge symmetries. Be- 
cause the higher spin symmetries are necessary for the consistency of the higher spin 
gauge theory, it is desirable to find a regularization which does not break them. Sec- 
ondly, because the heat kernel regularization depends explicitly on the background 
spacetime, it is not suitable for a nonperturbative study. For such a purpose, a 
regularization by finite- iV matrices, in which we can perform calculations directly in 
terms of matrices, is desirable. Thirdly, although the heat kernel regularization does 
not give the Einstein-Hilbert action, the Ricci-flat spacetimes are classical solutions, 
as we saw in the introduction. This is possible because in the derivation of the equa- 
tion of motion of the matrix model, Eq. (|1-12|) [or Eq. (|1-14|) ]. we varied not only 
the graviton but also all the fields with any spin. This can easily be done in terms 
of matrices, but once we rewrite the matrix model in terms of the local fields, such 
a calculation becomes hopelessly complicated; calculations become difficult already 
at the classical level. 



§3. Regularization by finite-iV matrices 



In the previous section we discussed the heat kernel regularization. Although 
it can be applied to any background spacetime and has an apparent physical in- 
terpretation, the actual calculation seems to be almost impossible, because in the 
regularized action, an infinite number of higher spin fields couple in a complicated 
manner. Furthermore, because it depends explicitly on the background spacetime, 
it is not suitable when we consider nonperturbative dynamics, e.g. the dynamical 
generation of the spacetime. These facts motivate us to introduce a regularization 
by finite-iV matrices. 
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In general, writing an explicit form of such a regularization is difficult, but there 
are classes of manifolds which have large symmetries, so that the corresponding 
Ep r i n have Lie group structures and allow explicit finite- iV regularizations that can 
be obtained using representation theory. In this section, we present a few examples: 
S n , RP n , CP™, T n and R n . We explain the case of S 2 in detail and show how 
the local fields and the generators of the diffeomorphism, the local Lorentz transfor- 
mation, and the higher spin gauge transformations are embedded in the matrices. 
Similar arguments can also be applied to other backgrounds. 

3.1. The covariant derivative on S n 

Let us consider the n-sphere S n with an isotropic and homogeneous metric. The 
commutation relations of the covariant derivatives V( a ) (a = 1, • • • ,n) are given by 

t v <°» v «i - {« J " «■'*'} °® - i£vf%- P' 1 ' 

where R is the scalar curvature. Here is the generator of Spin(n), which satisfies 
the commutation relation 

P { 2\0%] = \ {5 ac 0t$ - 5 bc O% - 6 ad o£ } + M®} • (3-2) 
It acts on V( a ) as 

Pt\ V (o) ] = \ {5 ab V (c) - 5 ac V (b) } . (3-3) 

Therefore, the algebra generated by V and is equivalent to spin(n + 1) under 
the identification 

V(a) ^ ^^ZTfinX Off — Ot +1 \ (34) 

where 0^ +1 ^ is the generator of Spin(n + 1). Therefore, E pr i n corresponding to S n 
is Spin{n + 1): 

E prm {S n ) = Spin{n + l). (3-5) 

This result is in some sense trivial, because it is well known that 

S n = Spin(n + 1) / 'Spin(n). (3-6) 

Here, Spin(n) is interpreted as the local Lorentz degrees of freedom. 

In order to express the covariant derivative in terms of a set of matrices, we must 
take the Hilbert space to consist of the smooth functions on E pr i n (S n ) = Spin(n+1), 
which is the regular representation, by definition. It is decomposed as 

C°° (Spin(n + 1)) = e r: i rr . rcpr . {V r • • • ® VX (3-7) 

dimv r times 
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where in the r.h.s., r runs over all the irreducible representations of Spin(n + 1). As 
we explain for the case of S 2 in detail below, by discarding the larger representations, 
we obtain a sensible cutoff. This is an analogue of the heat kernel regularization, 
because the Casimir operator of spin(n + 1) is simply the Laplacian on E pr i n (S n ) = 
Spin(n + 1). 

3.1.1. The covariant derivative on S 2 

The covariant derivative and a truncation of the Hilbert space 

As a concrete example, let us consider the case of n = 2 in detail. In this case, we 
can identify V and O with the generators Jj of Spin(3) = SU(2) as 

, 2 1 





<— ► 2 J3 ' (3-8) 
where generators Jj satisfy the commutation relations 

[Ji, Jj] = ie ijk Jk- (3-9) 
Let us label the states in the Hilbert space by the "angular momentum" (l,m) 



as 



S 2 \l, m; i) = 1(1 + l)\l, m; i), Js\l, m; i) = m\l, m; i). (3-10) 

Here, i = —I, — I + 1, • • • , I is the label for the different spin I representations in the 
regular representation*) . In our interpretation, the eigenvalue m of J3 is the spin of 

the state. This state has momentum yj (1(1 + 1) — m 2 ), because 

- V 2 \l,m;i) = I (f -Jl) \l,m;i) = - (1(1 + I) -m 2 ) \l,m;i). (3-11) 

For each spin and momentum, there are 21 + 1 states labeled by i . To understand 
the reason for this degeneracy, let us consider the case of a scalar field, m = 0. 
A scalar field on a sphere can be expanded in spherical harmonics. In terms of 
spherical harmonics, the state with momentum 1(1 + 1) is the state with the angular 
momentum /, which has 21 + 1 degrees of freedom. 

To see the correspondence discribed above more directly, it is instructive to 
recall the proof of the irreducible decomposition (|3-7|) . Let G be a compact Lie 
group. [ In the present case, we have G = Spin(3). ] We use the orthonormality of 
the representation matrices, expressed as 



\— J d 9 R$( 9 y4?(g) = j/ r){r,) S lk 5 jh (3-12) 



Vol(G) 

where dg is the Haar measure, (g) is the representation matrix for the irreducible 
representation r, and d r = dimV^- is the dimension of the representation r. By 
the Peter- Weyl theorem, the representation matrices R^'(g) form a complete set of 

*' We use i = — I, ■ ■ ■ ,1 instead of 1, • • ■ ,21 + 1, because it is related to the "angular momentum" , 
as we see below. 
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smooth functions from G to C. Therefore, any smooth function f(g) can be expanded 
as 

f(g) = Y.4 ) ^4 ) (g), (3-13) 

r 

where the sum is taken over all irreducible representations. We now see how the 

(r) 

quantities c\- transform under the action of G. The right-hand side transforms as 
h : $B$(g) -> $R${h-*g) = $ R% (h^)4] (g), (3-14) 

(r) 

which shows that c}/ transforms as 

l j 

$ - T 4 > (^ 1 ) C g- (3-15) 
The index i of Cjj transforms as the dual representation of r, while the index j is 

(r) 

invariant. It follows that c\- represents d r copies of the dual representation of r. 
Thus, we have found that the regular representation is decomposed as in Eq. (|3-7j) . 
In the present case, the state \l, m; i) corresponds to \/2Z + IK^Ag). Now consider 
the global right action 

m ■ f{g) f(gh), (3-i6) 

which commutes with the left action and acts on E pr i n (S 2 ) = Spin(3) as a rotation. 
[ Note that such a transformation exists only when E pr i n is a Lie group.] Under this 
transformation, the index i of \l,m;i) is rotated as follows: 

f(h) : \l,m;i) ^ Rf} (h)\l , m; j) . (3-17) 

Therefore, the states \l, 0; i) (i = —I, — I + 1, • • • , I) correspond to the spherical har- 
monics with angular momentum (I, i), while the states \l, m; i) are the corresponding 
modes of spin m fields. 

Next let us consider how to express the covariant derivative V( a ) in terms of a 
set of finite- N matrices. As explained above, the smaller representations correspond 
to the states with the smaller momenta and the lower spins. Therefore, it is natural 
to cut off the larger representations*) . If we discard the representations whose spins 
are larger than L, we have the representation 

v L = e m e • • • e v t ) 

1<L, Z:half integer v v ^ 



2/+1 times 

Span i |Z, m; i 



' = 0,^ ••• ,L; m,i = -l,--- ,1}. CMS) 



whose dimension is 



N = Yl (2/ + 1) 2 = -^(4L + 2)(4L + 3)(4L + 4). (3-19) 



24 

0<1<L, leZ/2 



*^ This can be regarded as a finite- N realization of the heat kernel regularization. 
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Because the momentum of the state \l, m; i) is y ^ W + 1) — wi 2 ), the momentum 
cutoff becomes more stringent as the spin m increases. The covariant derivatives can 
be expressed by the rotation generators acting on this space. In the following, we 
denote the projection to Vl by ttl, but we omit it for simplicity, unless it is confusing. 

Diffeomorphism, local Lorentz symmetry and local fields 

Because we are using a rather simple basis, diffeomorphisms and local Lorentz trans- 
formations can be written explicitly in terms of components. 

First, let us consider the operator /, which acts on the Hilbert space C°° (Spin(3)) 
as a multiplication by /(x, </):*) 

f:h(x,g)»f(x,g)-h(x,g). (3-20) 

Here, x and g denote the coordinates of M and Spin(2) = U(l), respectively. Such 
an operator can be written as a linear combination of Rmn in the form 

/=j;£.^, (3-2i) 

l,m,n 



where 



fid) = E f$n ■ V2lTlR<±(x,g). (3-22) 

l,m,n 



Because a product of representation matrices satisfies the relation 

*/2T+Tfl2„(x,s) ■ V2k + lR\f(x,g) 



k+l 



= Y, C(L;l,m,n;k,iJ)-V2l^R { X tn+j (x,g), (3-23) 

L=\k-l\ 

where C(L;l,m,n;k,i,j) is a product of Clebsch-Gordan coefficients, given by 
C(L;l,m,n;k,i,j) 

(2 ' + 2L + l +1) (<A^fc,m + ^,m)|fc,i)r(£,Z,fc,n + j|Z,n)|fc,j),(3-24) 

Rmn can be written as 

V21 + lR$ n = C(L;l,m,n-k,i,j)\L,m + i;n + j){k,i;j\. (3-25) 

L,k,i,j 

Because we are using the basis in which the covariant derivative and the Lorentz 
generator are block-diagonal, a multiplication of functions is basically off-diagonal. 



Although in the introduction we did not distinguish / and /, here we distinguish them so 
that / is not misinterpreted as an element of the Hilbert space. 



Regularization of the Covariant Derivative on Curved Space 11 

Note that R$ n raises the spin of the state by m, and hence it can be regarded as an 
operator of spin m.*' 

Next, let us consider the local Lorentz transformation. Because a parameter 
X(x) = A 01 (x) of the local Lorentz transformation is a scalar field on S 2 , it can be 
expanded in fl£>: 

X(x) = Y, Ai° • V5T+1r£(i), A = £ A» • ^2TTTi?£. (3-27) 

l,n l,n 

Using this, we can write the generator of the local Lorentz transformation as 

A = AJ 3 . (3-28) 

In the case of the diffeomorphism, a parameter X^(x,g) is a vector on S 2 , and 
therefore it can be expanded in R± x n : 



X^(x,g) = £ A™ ■ V2lTlR^ n (x,g), 

l.n 

AC* 1 ) = A™ ■ V2TTTRi{ n . (3-29) 



? (0 

y\ n ■ ■ ■ \ ^ -r xx 

l.n 

, n h ■ v ^ ~r xx 

l.n 

By using J± = J\ ± iJ2, the generator of the diffeomorphism becomes 

Similarly, we can express all the higher spin gauge transformations explicitly in terms 
of matrices. 

When we regularize the Hilbert space C°°(Spin(3)) using Vl, we must restrict 
A to act only on Vl- 

A — >A l = tt l Att l . (3-31) 

Note that this projection does not commute with the multiplication, 

{AA') L ^A L A' L , (3-32) 

but a pathology exists only near the momentum cutoff. 

The expansion in local fields given in Eq. (|1T1[) is expressed in terms of matri- 
ces, just as in the case of A. 

In Ref. 1), it is shown that the field f^(x,g) of representation r transforms under the right 
action f(h) : f(x,g) i— > f(x,gh) of Spin(n) as 

nh- 1 )fl]r(h)=RU (b) (h)f^. (3-26) 

In the present case, because Spin(2) — U(l) is abelian, the left and right actions are the same, and 
it is generated by the adjoint action of J3. Therefore, Rmn has spin m also in this definition. Note 
that the right action of Spin(2) considered here is different from that of Spin(3), defined by Eq. 
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Perturbative instability and its possible elimination 

Once the regularization is specified, we can calculate various quantities. Unfortu- 
nately, however, it is not possible to carry out a perturbative expansion about this 
background, because of the presence of the tachyonic modes, which come from the 
higher spin fields. This can be seen as follows. First, we decompose the matrix 
variables into the background and the fluctuations, as 

A a = P a + a a , ip a = + Xa , (3-33) 

where P a = iV( a ) = \p^J a for a = 1,2 and P a = otherwise, and add the gauge 

fixing and ghost terms Sgf+fp = —Tr (| [P a , a a ] 2 + [P a , b] [A a ,c]) . Then, the kinetic 
terms become 

S kin = -^Tr | ^a a (V 2 5 ab - 2iT ah ) a b - ^ X ?X + bV 2 c^ , (3-34) 

where T a = [P a , • ] and T ah = i[V a ,V b ]. In the present case, V 2 = f ( J? + J$) 
and T\2 = — ^Js, where Ji = [Jj, • ], commute and hence are simultaneously 
diagonalizable. The adjoint representation is labeled by the matrix indices (I, m; i) 
and (k,n;j), and each of them decomposes into the direct sum of (2k + 1)(2Z + 1) 
copies of spin \l — k\, ■ ■ ■ ,1 + k representations. These bases are related by the 
Clebsch-Gordan coefficients as 



X {l,m;i)(k,n;j) < > ^(s,p;i,j) = X (l,m;i)(k,n;j) x (h k, S, p\ 1, 171) \ k, n)' , (3-35) 

l,m;k,n 

where \k,n)' is the dual representation of \k,n). Here, Xi s ^a is a spin (s,p) state 



of J, which corresponds to a field of spin p *) and momentum y ^ (s(s + 1) — p 2 ). 

A straightforward but tedious calculation shows that adjoint representation we 
are considering can be decomposed as 



?s<2L 



(v s ®---®v s ), (3-37) 



k s times 

where 



| (s - 2L - 1) (2s 3 + 6s 2 + ALs 2 - 20sL -5s- 16L 2 s - 8L 2 - UL - 6) , (s = integer) 
^(2s + l)(2s + 8L + 7)(2s-4L-l)(2s-4L-3). (s = odd/2) 

(3-38) 



*' Strictly speaking, p represents the ^-dependence of at a ) bcd '" (x, g), which is obtained from 
ata)^ ® " ( x > d) in Eq. I1-1H by changing the indices with parentheses to those without: 

« W bcd '-X^,P)=^ (i '' ) (5)^ (c ' ) (3)^ ([i ' ) (p)---aK) (6 ' )<C ' )(ti ' ) '''^,3)- (3-36) 
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For simplicity, we have assumed here that L is an integer. Note that the dimension is 

Es=o( 2s + 1)** = (h( 4L + 2 )( 4L + 3 )( 4L + 4 )) 2 = ^ 2 > as expected. In this basis, 
we have 



{V 2 S ab - 2iT ab ) 



( s(s + l)-p 2 -2ip \ 

R 



2 



2ip s(s + 1) — p 2 



s(s + 1) — p 



2 



■ I 



ab 

(3-39) 



Smass = ~—Tr (A a A a ) , (3-41) 



but this is not positive-definite, because tachyons arise in the case 

(s(s + 1) - p 2 ) 2 - Ap 2 < 0. (3-40) 

Even if we introduce an infrared cutoff (that is, we discard small s), this inequality 
holds if p is sufficiently large. 

Some remarks are in order here. By adding a tachyonic mass term, 

R_ 

2ng 

we can make an n-sphere with scalar curvature R a classical solution *). There- 
fore, if such a term is generated dynamically, then the spacetime would compactify 
spontaneously. However, there remain some subtleties. First, the addition of such 
a term does not remove the tachyonic modes. Secondly, a fuzzy sphere and fuzzy 
torus also become classical solutions. The tree- level free energies on these back- 
grounds are of order — p-iV 3 , which are far smaller than those of ordinary spheres, 

~ — p-iV 5 / 3 . Therefore, the existence of tachyonic modes is not so surprising. If an 
ordinary sphere became stable in the large- N limit, then the dynamical correction 
would be very important: The higher spin fields should acquire mass dynamically 
and decouple from the low energy spectrum. It would interesting to study whether 
or not this is the case nonperturbatively, for example by employing a Monte Carlo 
simulation, as in Ref. 13). It would also be interesting to find other modifications of 
the action that make curved spaces classical solutions. Another way to remove the 
tachyons is to require A a to transform as a vector or scalar**) under the global right 
action of Spin(2): 

fih-^Aafih) = R a \h)A b {a = 1, 2), 

r(/i" 1 )^ a f(/i) = (a = 3, 4, • • • , 10). (3-42) 
Under this assumption, a a bcd "'(x,g), which is obtained from auy b ^ c ^^"' (x, g) in Eq. 

cm bv 

a a bcd -(x,g) = R a W(g)R b ( b '\g)R c ( c '\g)R d W(g) ■ ■ ■ a ia ,f^ d '>-(x,g),(3A3) 



In general, spaces satisfying the Einstein equation with cosmological constant ^j^-R become 
classical solutions. For details, see Appendix [B] 

**' Since a = 3, 4, • ■ ■ ,10 should be regarded as a label of scalar on S 2 , it is natural to require 
the second constraint in Eq. 13-42H . 
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does not depend on g. 1 ^ Because Spin(2) = U(l) is abelian, the left and right actions 
are the same, and it is generated by the adjoint action of J3. Therefore, imposing the 
above condition, the modes with \p\ > 1 are removed. Note that the modes removed 
here are those which do not allow simple interpretations in terms of low energy field 
theory. Although in the present case the constraint Q3-42j) can easily be written in 
terms of matrices as 

[J 3 , Ah] = ±A ±X (A± = A ± iA 2 ), [J 3 , A a ] = (a = 3, 4, • • • , 10), (3-44) 

on gereric backgrounds it is difficult to write this condition explicitly. Note that 
such a condition depends explicitly on the background. In any case, if we remove 
tachyonic modes by modifing the action or imposing the condition ()3-42j) . we can 
perform a perturbative calculation around it. As we saw above, the calculation in 
this case is almost parallel to that on a fuzzy sphere. As an example, we calculate 
the one-loop free energy under the constraint (|3-42j) . 

One-loop free energy under the constraint (|3-44|) 

We impose on ^,6 and c conditions similar to the constraint (|3-44j) . 

[J 3 , V] = ±^7 12 V>, [J3, b] = [J3, c] = 0. (3-45) 

Then, from Eq. (|3-34j) . the one-loop free energy becomes 

Fone-ioo P = \Tr b log (V 2 5 ab - 2iT ah ) - ^Tr f \ogf - Tr gh V 2 , (3-46) 

where the subscripts on Tr indicate that we take the trace only over the states 
satisfying Eqs. (j3-44|) and ()3-45[) . Each term in Eq. ()3-46|) is evaluated as follows: 

l -Tr h log (V 2 5 ab - 2iF ab ) 

= 5 ^ fc s log-s(s + l) + - 2^ k s log> 1 



2 ^ I s 2 (s + l)< 

simteger s:mteger ^ ' 



(3-47) 



~Tr f logf = --Tr f \og { V 1 -\ - -f 




= -2 £ k s log^\(s(s + l)-\) (3-48) 

s=odd/2 

-Tr gh V 2 = - Y, fc»log|*(a + l). (3-49) 

s:integer 

Therefore, we find 



Fone-loop = 4 Y ks lo S s ( s + !) - 2 ks log \ [ s ( s + X ) 

s:integer s=odd/2 



I r _ 1 

4 



Regularization of the Covariant Derivative on Curved Space 



15 



1 f (s(s + l) — l) 2 — 4 1 

+ 2 £ M°g s2 J + J 2 ■ (3.50) 

s: integer ^ J 

In order to remove tachyonic modes completely, we must introduce an infrared cutoff; 

that is, we should discard the modes satisfying s < 1 by hand. Then, numerically 
we find that it behaves as 

Fane-loop ~ ~L* ~ — A 4 / 3 , (3-51) 

and is negligible compared with the tree-level free energy. 



Difference from the fuzzy sphere 

Here we stress again that although V is expressed in terms of the generators of 
SU(2), it is completely different from the fuzzy sphere. We can express the ordinary 
two-sphere using two matrices, while in the case of the fuzzy two-sphere, three ma- 
trices are needed. Because the covariant derivative on the ordinary sphere is only 
a part of su(2), it is not a classical solution of the matrix model with the cubic 
term, 14 ) although the fuzzy sphere is a classical solution. ( For details, see Appendix 
IbI ) The radii are also different; the radii of the commutative and fuzzy spheres are 
~ -j= and -4=, respectively. 

3.2. The covariant derivative on MP' 1 

MP™ is obtained from S n by identifying antipodal points. The corresponding 
principal bundle E pr i n (MP n ) is also obtained through this identification, and so we 
have 

Epr in (RP n ) = Spin{n + l)/± = SO{n + 1). (3-52) 
This can be seen also from the relation 



SO{n + 1)/Spin(n) = RP n . (3-53) 

As in the case of Spin{n + 1), C°°(SO(n + 1)) can be expanded in the irreducible 
representations of SO(n + 1), and the restriction to smaller representations provides 
a natural cutoff, which is an analogue of the heat kernel regularization. 

3.2.1. The covariant derivative on RP 2 

Functions on SO(3) can be expanded in the representation matrices with integer 
spins: 

C°° (SO{3)) = © (Vi © • • • © Vi) 

l : integer s v ' 



22+1 times 



Span I \l, to; i) 



I = 0, 1, 2, • • ■ ; m, i = — I, — I + 1, 



Span | Rmi® I : integer; m,i = —I, —/ + !,••• , Z j . (3-54) 
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By restricting the value of I, we can introduce an explicit finite- N regularization, 



V L = © (Vi © • • • © Vi) = Span \ \l,m;i) 

1<L, I : integer v v ' \ 

2Z+1 times 



where 



I = 0, ,L; m,i = -/,••• ,1 } 



(3-55) 



N = -(2L + l)(2L + 2)(2L + 3). (3-56) 
6 

The covariant derivative is expressed in SO(3) generators acting on this space, just 
as in Eq. ()3-8|) , The only difference from the case of S 2 is that in the present case, 
only the integer spins appear*). Therefore, the local fields and various symmetries 
are embedded in the same way, and if we expand the matrices about this background, 
tachyonic modes appear, due to the existence of higher spin fields, just as in the case 
of S 2 . 

3.3. The covariant derivative on CP™ 

In the case of CP™, the covariant derivative and the Lorentz generators form 
a subalgebra of spin(2n + 2), and a procedure similar to that explained above can 
be applied. Although E pr i n (CP n ) is not so simple, if we are interested only in 
holomorphic quantities then the covariant derivatives on CP™ takes simple forms. 
Using the indices J and J, defined by 

V J = yV- 1 + y/=lv 2j , V J = V 2j ^ - V^lv 2j , (3-57) 

under the identification 



( > ]j n(n + l) J>+i ' lJj V n(n + l) W +M U u 

(3-58) 



and 



n+l 



T£off +2) =0, (3-59) 

7=1 



the algebra generated by V becomes su(n + 1). Therefore, taking C°° (SU(n + 1)) 
as a Hilbert space, the covariant derivative can be expressed in terms of the rotation 
generators acting on it. C°° (SU(n + 1)) can be decomposed into the irreducible 
representations of SU(n + 1), and regularization by finite- iV matrices is realized by 
cutting off the larger representations. 



*' Rij^ is symmetric (resp., antisymmetric) under the identification of the antipodal points of 
Spin(3) ~ S" 3 if / £ Z> (resp., I G Z> + ±). Functions on RP 3 = 50(3) are the functions on 
S 3 = Spin(3) that are symmetric under this identification. 
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3.4. The covariant derivative on IR n and T n 

Next we consider a flat space IR n and a torus T n . In this case, no tachyon 
appears, and a perturbative calculation is possible without any constraint. 

The covariant derivative and the Lorentz generators are identified with the trans- 
lation and rotation generators of the Poincare group ISO(n). Therefore, the princi- 
pal Spin{n) bundle corresponding to W 1 is the Poincare group: 

E prm (R n ) = ISO(n). (3-60) 

Compactifying E prin (M. n ), we obtain E prin (T n ). 

Although EpriniW 1 ) has a group structure, for actual calculations it is better to 
use another regularization. The strategy is to use a momentum basis in the direction 
of W a (or T n ), while in the direction of Spin(n) we use a coordinate basis*). 15 ) With 
such a basis, the covariant derivative can be written as 

(V(o)) jj = iP(a) S V = i (^(a) 6 ^ 1 )) (PbT 5af30~ij, (3-61) 

where i,j and a, (5 indicate the degrees of freedom of W 1 and Spin(n), respectively, 
and we have / = (i,a) and J = (j,(3). The quantity p\ is an eigenvalue of — idb, 
which is the momentum in the direction of M. n . We assume that p\ is uniformly 
distributed below the momentum cutoff. Then, p 1 ,, is also uniformly distributed. 
Therefore, by regarding p^ as the "momentum along the (a)-th direction" , we can 
apply the techniques of the quenched reduced model. 16 ) 

§4. Conclusions 

In this paper, we discussed a method for regularizing the covariant derivative 
on a Riemannian manifold by a set of finite- iV matrices. Such explicit constructions 
are important for various reasons. As an example, consider spontaneous spacetime 
generation. 17 ) In our interpretation, the four-dimensional eigenvalue distribution 
observed in Ref. 18) suggests that the four-dimensional spacetime is realized dy- 
namically. However, extensity of the eigenvalues does not reveal the topology **). 
(Note that any manifold entails infinitely large momentum.) In order to extract 
information concerning the topology, we must study the detailed structure, such as 
commutation relations and the degeneracy of eigenvalues. 

Although for generic manifolds it is difficult to write down explicit forms of 
regularizations, for some classes of manifolds that have higher symmetries, we can 
introduce explicit cutoffs, and explicit calculations are possible. We studied spheres, 
real and complex projective spaces, tori and flat spaces as examples. 

In the case of the n-dimensional sphere S n and the real projective space ILP n , the 
corresponding principal Spin(n) bundles are Spin(n+1) and SO(n+l), respectively, 
and the covariant derivative can be expressed in terms of the rotation generators 
acting on the regular representation. By discarding the larger representations, we 

In such a basis, fields with various spins are mixed. 
**) This is in a sharp contrast to the original interpretation. 17 ' 
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can introduce a cutoff for both the momentum and spin of the fields. The meaning 
of the states in the Hilbert space is clear. The covariant derivatives on S n and MP 11 
satisfy the same commutation relation, but they are expressed in terms of different 
matrices; indeed, the degeneracies of the irreducible representations are different. 
This is a reflection of the fact that S n and MP n are locally isomorphic but globally 
non- isomorphic. 

We investigated the case of S 2 in detail and found how the local fields and the 
generators of the diffeomorphism, the local Lorentz transformation, and the higher 
spin gauge transformations are embedded in the components of the matrices. If we 
choose a Yang-Mills-type action, such as that of IIB matrix model, a perturbative 
expansion about this background cannot be carried out, due to the existence of 
tachyonic modes. Such modes come from the local fields, which depend on the 
coordinate of Spin{2). If we use an action that does not entail such tachyons, then 
using the techniques presented in the calculation becomes almost parallel to that 
on a fuzzy sphere background. 

In the case of a flat background, we can use the techniques of the quenched 
reduced models. 16 ) In this case, no tachyonic modes appear. Because the quenched 
reduced model reproduces the U(N) gauge theory, by studying this background, we 
may find new connections between gauge theory and gravity. A calculation of the 
correlators of the vertex operators 20 ) would be very useful in this context. It would 
also be useful to clarify the relation between our formalism and that for the original 
interpretation of IIB matrix model. 
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Appendix A 

The Heat Kernel Regularization 



In this appendix, we present the formulae used in £J2I We start with the 
well-known formula 

Tr t l = n d ~ D J dg J ed d x(x, g\e tAprin \x, g) 

d-D t -D/2 r r / t \ 
= {An)D/2 J dgj ed d x (l + -R prm + 0(t 2 )J , (A-l) 

where R pr i n is the scalar curvature of E pr i n , and D = i s the dimension of 

E pr i n . At low energies, we can ignore all terms but the first in the limit t — ► 0. 

For simplicity, let us denote both the subscript (a) of V( a ) and ab of O a b by a 
single index, I, and denote V( a ) and nO a b by D\. Then, we can evaluate the trace 
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of f IJ "' K (x,g)DiDj ■ ■ ■ Dk (we assume I,J, -- , K to be symmetric) as follows. 
By the general covariance, the trace of f IJ (x,g)DjDj takes the form 

Tr t (f IJ (x,g)D I D J ) 

= K d - D t' D l 2 ' x f dg f ed d x (a5uf IJ + ptR prm f IJ 5u + jt(R pnn )uf IJ + •••), 

(A-2) 

where R W in and (R pr in)u are the scalar curvature and the Ricci tensor of E pr i n (M), 
respectively. The coefficients can be determined by taking f IJ = 5 IJ . For example, 
from the relations 



Tr t {S IJ DjDj) = Tr t {\ rin ) 

= Jt^^ 

D K d - D t- D l 2 - 1 



2 (4vr)^/ 2 



Q dg J ed d x^j +■■■ (A-3) 



and 



Tr t {5 IJ DiDj) = K d ' D r D l 2 - x J dg J ed d xa5 IJ 5 IJ + ■■■ 



we have 



a = ~2^W (A ' 5) 

Similarly, we have 

Tr t {f IJKL (x,g)D I D J D K D L ) 

= K \^J/ 2 2 j d 9 j ed d xf IJKL [SijSkl + SikSjl + 5 IL S JK ) + • • • , (A-6) 
Tr t (f IJKL (x,g)D I D J D K D L D M D N ) 

K d-D t -D/2-3 r r 

= 8 ^D/2 J dg J eddx f IJKL (5u s klo~mn + 14 permutations) + • • • . 

(A-7) 

Note that terms with an odd number of Dj vanish. 

Appendix B 

Deformed Actions and Their Classical Solutions 



As we saw in the Introduction, with the ansatz 

A a = iV (a) , Va = 0, (B-1) 
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the equation of motion of IIB matrix model gives Ricci flat spacetimes as classical 
solutions. Note that Ricci-flat spacetimes with fewer than 10 dimensions are also 
classical solutions; indeed, the only difference is that in this case we use the ansatz 

A a = iV {a) {a = !,-■■ ,n), A a = (a = n + 1, • ■ • , 10), ^ a = 0. (B-2) 

By deforming the action, many noncommutative spaces become classical solu- 
tions. In this appendix we determine the kinds of classical solutions that exist in our 
interpretation. 

B.l. Matrix model with a mass term 
If we add a mass term 



m 

-—Tr(A a A a ) 

g 2 



to the equation of motion becomes 



A h 



An. -A 



+ 2m 2 A a = 0. 



(B-3) 



(B-4) 



Here we set ip = for simplicity. Substituting the ansatz (|B-1|) or (|B-2|) into (|B-4|) 
we have 



V a R ab cd = 0, R ab = 2m 2 5 ab . 



(B-5) 



Using the Bianchi identity, the former follows from the latter. Note that Eq. (jB 
is the Einstein equation with a cosmological constant 

(n — 2)m 2 . 



A 

B.2. Matrix model with a cubic term 

Next, let us consider the matrix model with the cubic term 14 ) 



S, 



cs 



3<? 2 



e abc Tr([A a ,A b ]A c ), 



(B-6) 



(B-7) 



which has the fuzzy sphere as a classical solution. Here e abc is totally antisymmetric, 
g i23 _ ]_ an( j e abc = o if at least one of the indices is not 1,2 or 3. 



The equation of motion is 
~A b 



A a , A 



iae abc [A b ,A c ] = 0, 



and using the ansatz (|B-1(I or ()B-2|) . we obtain 



Rab Cd = 



(B-8) 

(B-9) 
(B40) 
(B41) 

In this way, it is found that with the ansatz (jB-lj) and ()B-2j) . classical solutions of 
the matrix model with the Chern-Simons term are rather trivial. 



for n = 2, 3 and 

V b R ab de - ae abc R bc de = 0, R ab = 
for n > 4. Using the Bianchi identity, we can rewrite the latter as 



Rab — U, e rc bc — U. 
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